function rest=testfunc(t,y)



% rest=zeros(termM,1);
rest(1,1)=-y(2);
rest(2,1)=y(1);
rest(3,1)=y(3);
% rest=[y(2);y(1)];
end
